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Abstract.  The  activity  patterns  of  rhythmically  firing  neurons  in  monkey  primary  somatosensory 
cortex  (SI)  were  studied  during  trained  wrist  movements  that  were  performed  , 

palmar  vitoation.  Of  1,222  neurons  extracellularly  recorded  in  SI,  129  cells  (-11%)  discharged 
rhvthmically  (at  -30  Hz)  during  maintained  wrist  position.  During  the  initiation  of  vibratory-cued 
movements,  neuronal  activity  usually  decreased  at  -25  ms  after  vibration  onset  followed  by  an 
additional  decrease  in  activity  at  -60  ms  prior  to  movement  onset,  ^ythmic^ly  firing  neurons  are 
not  likely  to  be  integrate-and-fire  neurons  because,  during  activity  changes,  their  rhythmic  firing 
pattern  was  disrupted  rather  than  modulated.  The  activity  pattern  of  rhythmwally  firing  neurons 
was  complimentary  to  that  of  quickly  adapting  SI  neurons  recorded  during  the  perfoimance  of  this 
task  INelson  et  ah,  1991).  Moreover,  disruptions  of  rhythmic  activity  of  individual  SI  neurons 
were  similar  to  those  reported  previously  for  local  field  potential  (LFP)  oscillations  m  sensonmotor 
cortex  during  trained  movements  (Sanes  and  Donoghue,  1993).  However  rhythmic  activity  of  SI 
neurons  did  not  wax  and  wane  like  LFP  oscillations  (Murthy  and  Fetz,  1992;  Sanes  and 
Donoghue,  1993).  It  has  been  suggested  that  fast  (20  -  50  Hz)  cortical  oscillations  he  initiated 
by  inhibitory  intemeurons  (Cowan  and  Wilson,  1994;  Llinas  et  al.,  1991;  Stem  and  Wilson, 

1994).  We  suggest  that  rhythmically  firing  neurons  may  tonically  inhibit  quickly  adapting  neurons 
and  release  them  from  the  inhibition  at  go-cue  onsets  and  prior  to  voluntary  mov^ents.  It  is 
possible  that  rhythmically  active  neurons  may  evoke  interrmttent  oscillations  in  other  cortical 
neurons  and  thus  regulate  cortical  population  oscillations. 

Introduction 

During  the  execution  of  motor  tasks,  large  groups  of  sensorimotor  cortical  neurons  mtermittently 
become  involved  in  coherent  rhythmic  activity  (Murthy  and  Fetz,  1992;  Sanes  and  Donoghue, 

1993)  Rhythmic  activity  has  been  suggested  to  be  a  functionally  important  type  of  neuronal  tiring 
that  may  serve  to  switch  between  behavioral  modes  and  to  establish  dynamic  coupling  between 
cortical  areas  (for  review,  see  Gray,  1994;  Llinas,  1990;  Lopes  da  Silva,  1991;  Sheer,  9  , 

Singer,  1993;  Steriade,  1993).  Rhythmic  activity  in  the  form  of  local  field  potential  (LFP) 
oscillations  at  -40  Hz  has  been  studied  most  extensively  in  cat  visual  cort^  (^khorn  et  al.,  19  , 

Gray  and  Singer,  1989;  for  review,  see  Engel  et  al.,  1992;  Gray  et  al.,  1991;  Gray,  19^;  Singer, 
1993).  Similar  LFP  oscillations  have  been  demonstrated  in  sensorimotor  cortex  in  cat  (Bouyer  et 
al.  1981  1987)  and  in  monkey  (Murthy  and  Fetz,  1992;  Rougeul  et  al.,  1979;  Sanes  and 
Donoghue,  1993).  In  these  studies,  it  was  suggested  that  oscillations  in  motor  and  sensory  cortical 
areas  play  a  role  during  behavior  in  sensorimotor  integration  and  focal  attention. 
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Individual  neurons  having  oscillatory  properties  may  be  important  for  the  initiation  of 
cortical  oscillations  (Llinas,  1990;  Steriade,  1993).  Single  morphologically  identified  neurons  with 
40-Hz  rhythmicity,  probably  inhibitory  intemeurons,  have  been  recorded  in  cortical  slice 
preparations  (Llinas  et  al.,  1991).  The  presence  of  fast  (~40  Hz)  rhythmic  IPSPs  has  been 
demonstrated  while  recording  from  pyramidal  cortical  neurons  in  vivo  in  urethane-anesthetized 
rats  (Cowan  and  Wilson,  1994;  Stem  and  Wilson,  1994).  These  rhythmic  IPSPs  may  have 
resulted  from  cortical  intemeuronal  inputs.  In  addition,  depolarization-induced  20-  to  40-Hz 
oscillations  have  been  reported  for  a  subset  of  cortical  neurons  with  long  axons  recorded  in  vivo 
in  urethane-anesthetized  cats  (Nunez  et  al.,  1992).  Therefore,  both  intrinsic  properties  and 
extrinsic  circuitry  may  be  involved  in  generating  rhythmic  activity. 

While  the  majority  of  data  demonstrating  oscillations  in  sensorimotor  cortical  areas  has 
been  obtained  while  recording  LFPs  during  behavior  (Bouyer  et  al.,  1981;  1987;  Murthy  and  Fetz, 
1992;  Rougeul  et  al.,  1979;  Sanes  and  Donoghue,  1993),  some  studies  have  examined  rhythmic 
activity  of  single  neurons.  Murthy  and  Fetz  (1994)  described  monkey  sensorimotor  cortical 
neurons  that  intermittently  generated  rhythmic  discharges  during  episodes  of  LFP  oscillations. 
Ahissar  and  Vaadia  (1990)  reported  that  second  somatosensory  cortical  (SII)  neurons  recorded  in 
awake  monkeys  often  exhibited  rhythmic  activity  that  was  sustained  during  immobility.  The 
pattern  of  rhythmic  firing  of  these  SII  neurons  was  disrupted  by  tactile  stimulation  or  during 
voluntary  movements.  Ahissar  and  Vaadia  (1990)  suggested  that  oscillatory  activity  of  SII 
neurons  may  be  important  for  texture  analysis.  The  functional  role  of  rhythmically  active  neurons 
in  the  primary  somatosensory  cortex  (SI)  of  behaving  animals  is  not  completely  understood. 

To  determine  the  changes  in  activity  patterns  of  rhythmically  firing  SI  neurons  during  the 
initiation  of  somatosensory-cued  movements,  we  analyzed  data  from  experiments  in  which 
monkeys  performed  voluntary  wrist  flexions  and  extensions  in  response  to  vibrotactile  go-cues. 
We  examined  the  characteristics  of  rhythmic  spike  trains  during  a  period  when  vibration  was  not 
present  as  well  as  how  these  spike  trains  changed  in  association  with  the  presentation  of 
vibrotactile  go-cues  and  movement  initiation.  Also,  we  sought  to  relate  the  activity  of  rhythmically 
firing  neurons  with  that  of  other  SI  neurons  recorded  during  the  same  experimental  paradigm 
(Nelson,  1988;  Nelson  et  al.,  1991;  Lebedev  et  al.,  1994).  Based  on  these  observations,  we 
constructed  a  model  depicting  how  rhythmically  firing  neurons  may  influence  the  firing  patterns  of 
other  cortical  neurons.  Some  of  the  data  have  been  presented  previously  in  preliminary  form 
(Lebedev  and  Nelson,  1993). 

Methods 

Experimental  Apparatus  and  Behavioral  Paradigm 

Six  adult  male  rhesus  monkeys  (Macaca  mulatta;  monkeys  C,  F,  G,  H,  M,  and  N)  were  trained  to 
perform  sensory-triggered  wrist  movements.  The  monkeys  were  cared  for  in  accordance  with  the 
NIH  Guide  for  Care  and  Use  of  Laboratory  Animals,  revised  1985.  Each  animal  sat  in  an  acrylic 
monkey  chair  with  its  right  forearm  on  an  armrest  and  its  right  palm  on  a  moveable  aluminum  plate 
(Fig.  lA).  One  end  of  the  plate  was  attached  to  the  axle  of  a  brushless  DC  torque  motor  (Colburn 
and  Evarts,  1978).  A  load  of  0.07  Nm  was  applied  to  the  plate,  which  assisted  wrist  extensions 
and  opposed  flexions.  Thus,  to  maintain  a  centered  wrist  position,  the  monkeys  had  to  rnaintain 
pressure  upon  the  plate  using  the  palmar  surface  of  their  hands.  Feedback  of  wrist  position  was 
provided  by  a  visual  display  located  35  cm  in  front  of  the  animals.  This  display  consisted  of  3 1 
light-emitting  diodes  (LEDs)  (Fig.  IB).  Current  wrist  position  was  indicated  by  illurmnating  one 
of  the  LEDs.  The  middle,  red,  LED  corresponded  to  a  centered  wrist  position.  Yellow  LEDs 
above  and  below  the  middle  LED  indicated  successive  angular  deviations  of  1  °. 

Each  monkey  was  trained  to  make  untargeted  ballistic  wrist  flexions  or  extensions  in 
response  to  vibrotactile  stimulation  of  their  palms  through  the  plate  (Nelson,  1988;  Nelson  et  al., 
1991).  Movements  were  executed  by  the  same  hand  that  was  stimulated.  The  experimental 
paradigm  is  illustrated  in  Fig.  1C.  To  begin  a  trial,  the  monkey  first  centered  the  plate.  At  this 
time,  a  movement  direction  request  was  given  by  the  presence  or  absence  of  illumination  of  an 
instructional  red  LED.  This  LED  was  located  in  the  upper  left  corner  of  the  visual  display  (Fig. 
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IB).  Illumination  of  this  LED  signaled  that  extension  was  the  appropriate  movement.  If  this  LED 
was  not  illuminated,  the  appropriate  movement  was  flexion.  The  monkey  was  required  to  hold  the 
SeTn  theTnteredpositiLforO.5, 1.0,  1.5,  or  2.0  s  (pseudorandomzed).  Movements  of  mom 
than  0  5°  from  the  center  during  this  hold  period  canceled  the  trial.  After  the  monkey  pccessfully 
completed  the  hold  period,  the  plate  was  vibrated  by  driving  the  torque  motor  with  a  sine  wave  at 
27  57  or  127  Hz  The  angular  deflection  of  the  handle  during  the  vibration  was  less  than  0^6  . 
The  onset  of  vibratory  stimulation  thus  served  as  a  cue  for  movement.  When  a  movement  of  at 
least  5°  in  the  required  direction  was  made,  the  vibration  was  turned  off,  and  the  animal  received  a 
fruit  juice  reward.  A  new  trial  began  when  the  animal  once  again  centered  the  plate. 

Electrophysiological  Recordings  and  Histology 

Once  an  animal  achieved  stable  performance  in  the  task  (-2,000  rewarded  trials  per  experiment^ 
Ssfonra  stainless  steel  recording  charnber  was  surgically  implanted 
extracellular  recordings  of  the  activity  of  SI  neurons  (see  Nelson  et  al.,  1991,  for  details). 
Transdural  penetrations  were  made  daily  into  the  region  of  SI  that  represented  the  contralateral 
hand.  Platimm-iridium  microelectrodes  with  impedances  of  1  -  2  MOhms  were  used  for 
recordings.  The  depth  of  the  electrode  was  varied  to  achieve  the  best  isolation  between  the 
neuronafwaveform  in  question  and  other  visible  neuronal  activity.  The  activity  of  single  units  was 
amolified  and  discriminated  using  a  time  window  discriminator  having  two  independently 
controllable  window  thresholds.  Each  channel  (upper  and  lower  window)  of  the  discrinunator  w^ 
armed  by  a  waveform’s  first  crossing  of  that  threshold.  An  acceptance  pulse  of  0.2  ms  was  issued 
when  tSt  threshold  was  recrossed.  The  time  from  initial  threshold  crossing  to  the  end  of  the 
acceptance  pulse  was  typically  on  the  order  of  0.5  ms.  As  well,  the  electronic 
the  entire  data  collection  system  while  running  the  paradigm  was  confimed  to  be  on  the  order  of 
0.4  ms  by  analyzing  recordings  of  white  noise.  Pulse  data  was  stored  in  a  computer  by 

conventional  means  (Evarts,  1966;  Lemon,  1984).  tacV  hw  HoViHv 

Neuronal  receptive  fields  (RFs)  were  manually  exaimned  outside  the  task  by  lightly 
touching  skin  surfaces  with  hand-held  probes,  manipulating  joints  and  palpating  muscles^  An  1^ 
was  classified  as  “cutaneous”  if  the  neuron  preferentially  responded  to  light  touch,  and  as  deep  if 
the  neuron  responded  to  bending  a  joint  or  to  muscle  palpation.  To  ensure  that  stimulation 
overlying  skin  was  not  mistaken  as  a  response  from  a  muscle,  the  skin  was  displaced  laterally  and 

the  muscle  again  palpated.  For  some  neurons,  no  clear  RE  could  be  found. 

At  regular  intervals,  the  EMG  activity  of  forearm  muscles  acting  across  the  wrist  was 
recorded  (Nelson,  1987).  Intramuscular  EMG  wires  (stranded  stainless  steel,  TEFLONO 
insulated-  Bergen  Wire  Rope  Co.)  were  temporarily  implanted  in  muscles  by  using  sterile  25-gauge 
needles  as  guides.  EMG  activity  was  converted  into  pulse  data  with  a  window  discriimnator 
rVaadia  et  al.,  1988;  Lebedev  et  al.,  1 994)  and  stored  in  the  same  form  as  the  neuronal  data.  . 
^  On  the  last  recording  day,  electrolytic  lesions  ( 10  [lA  of  current  for  10-20  s)  were  made  in  ■ 
the  cortex  to  mark  locations  of  interest.  The  animals  were  then  deeply  anesthetized  with  sodium 
pentobarbital  and  transcardially  perfused  with  10%  buffered  formol-salme.  Histological  sagittal 

sections  of  the  cortex  were  prepared,  and  recording  sites  were  reconstructed  based  upon  the  depth  . 
,  ,  1  ^  _ 1  trk  tVip*  mjirlciTia  i  rl?.  2)  (JNClSOn,  ' 
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1988;  Nelson  et  al.,  1991). 


Analyses  of  Rhythmic  Spike  Trains  (  1 

Rhythmically  firing  SI  neurons  were  selected  by  examining  the  500-ms  epoch  of  the  spike  trmns 
that  immediiely  preceded  vibration  onset.  This  corresponded  to  the  shortest  duration  of  the  hold 
period  during  which  an  animal  actively  maintained  a  stable  wrist  position.  Because  of  this  time 
limit,  slow  rhythmic  activity  (<  8  Hz)  could  not  be  detected.  Conclusions  on  whether  neurons 
fired  rhythmically  or  nonrhythmically  were  made  by  exaimning  the  expectation  density  (ED) 
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histograms  (Mountcastle  et  al.,  1969,  1990).  ED,  also  termed  the  “autocorrelation  function” 
(Zadeh,  1957),  represents  the  probability  of  a  neuronal  discharge  at  a  certain  time  after  a  given 
discharge.  For  rhythmic  spike  trains,  ED  histograms  contain  peaks  at  the  multiples  of  the  rhythmic 
period  (Fig.  3A)  (Perkel  et  al.,  1967;  Poggio  and  Viernstein,  1964;  Segundo  et  al.,  1968).  The 
average  height  of  subsequent  bins  in  the  ED  histogram  gradually  flattens  to  a  stabilized  level  (Fig. 
3A).  We  quantified  the  degree  of  rhythmicity  as  the  ratio  of  the  height  of  the  first  ED  peak  to  the 
stabilized  level  of  ED  histogram: 

E']=Ej/R,  (1) 

where  E'j  is  the  normalized  height  of  the  first  peak  measured  from  zero  level,  Ej  is  the  height  of 
the  first  peak,  and  R  is  the  stabilized  level.  The  stabilized  level  was  calculated  by  averaging  ED 
bins  in  the  range  0  to  200  ms. 

The  selection  of  rhythmically  firing  neurons  was  made  using  four  criteria:  (1)  mean  firing 
rate  of  the  selected  neurons  during  the  hold  phase  of  the  paradigm  exceeded  10  spikes/s;  (2  and  3) 
the  heights  of  the  first  and  second  ED  peaks  were  greater  than  1.5/?  and  1.2/? ,  respectively  (Fig. 
3C,  D);  (4)  time  intervals  from  zero  to  the  first  peak  {tj.  Fig.  3A)  and  from  the  first  to  the  second 
peak  {t2  -  tj)  differed  by  less  than  15%  (percentage  Jitter,  Fig.  3D).  Criterion  1  was  used  to 
initially  select  the  population  of  tonically  active  neurons.  Criterion  2  was  chosen  because  the 
distribution  of  E'j  for  tonically  active  neurons  indicated  a  major  subpopulation  with  in  the 
range  of  nonrhythmic  firing  (1  -  -1.5)  and  a  smaller  subpopulation  of  the  neurons  with  >  1.5 
(Fig.  3B  -  D).  Criteria  3  and  4  previously  have  been  used  by  Ahissar  and  Vaadia  (1990)  and 
confirmed  the  presence  of  multiple  peaks  in  ED  histogram.  The  distributions  of  E'j  for  the 
selected  groups  of  rhythmically  and  nonrhythmically  firing  neurons  were  statistically  different  (p 
<  0.0001;  unpaired  Mest),  and  had  a  small  overlap  in  the  range  of  E'j  from  1.5  to  -2.0  (Fig. 

3B).  For  rhythmically  active  neurons,  their  rhythmic  frequency  was  calculated  as 

/  =  2  /  [^;  +  (^2  -  ^/)])  (2) 
where  f  is  rhythmic  frequency  and  tj,  t2  are  the  intervals  from  zero  to  the  first  and  second 
peaks,  respectively.  Some  ED  histograms  showed  qualitative  features  that  were  characteristic  of 
bursty  activity  pattern  (Wilson  et  al.,  1977;  Wilson  and  Groves,  1981).  These  ED  histograms 
indicated  high  probability  of  firing  at  very  short  intervals  after  spike  generation  (<10  ms).  Firing 
probability  then  decayed  gradually  to  the  stabilized  level  over  a  period  of  20  -  50  ms.  Bursty 
neurons  were  excluded  from  these  analyses. 

Renewal  density  (RD)  histograms  and  joint  interval  plots  were  calculated  to  determine  the 
degree  of  serial  dependency  in  the  spike  trains.  We  sought  to  determine  whether  the  rhythmic 
spike  sequences  of  SI  neurons  were  reset  after  the  occurrence  of  each  discharge  (that  is,  could  be 
described  as  a  renewal  process)  (Perkel  et  al.,  1967)  or  more  complex  serial  dependencies  of  the 
ISIs  were  present.  Nonrenewal  spike  trains  often  occur  when  a  neuron  is  rhythmically  driven  by 
external  inputs  (Mountcastle  et  al.,  1969,  1990;  Surmeier  and  Towe,  1987a,  1987b).  Joint  interval 
scattergrams  display  each  ISI  as  a  function  of  its  immediately  preceding  ISI  (Fig.  4B).  Joint 
interval  methods  have  been  used  to  analyze  serial  dependencies  of  the  ISIs  (Rodieck  et  al.,  1962; 
Surmeier  and  Towe,  1987a,  1987b;  Siebler  et  al.,  1991).  Another  technique  for  detecting  serial 
dependencies  of  the  ISIs  is  the  RD  method  (Mountcastle  et  al.,  1969,  1990;  Perkel  et  al.,  1967; 
Poggio  and  Viernstein,  1964).  RD  is  the  ED  calculated  for  randomly  shuffled  ISI  sequences  (Fig. 
4D,  E).  For  a  simple  renewal  process,  ED  and  RD  histograms  are  not  different.  However,  if 
serial  dependencies  of  the  ISIs  are  present,  these  histograms  may  be  different.  Thus,  in  the  case  of 
externally  driven  activity,  a  peak  in  a  RD  histogram  is  smaller  than  the  corresponding  peak  in  the 
ED  histogram  (Mountcastle  et  al.,  1969,  1990).  To  estimate  the  difference  between  ED  and  RD 
histograms,  we  compared  the  magnitudes  of  the  first  and  second  peaks  in  the  ED  histogram 
normalized  by  the  stabilized  level  with  those  of  corresponding  peaks  in  the  RD  histogram. 

Analyses  of  Changes  in  Neuronal  Activity  and  EMG 

Changes  in  neuronal  activity  associated  with  go-cues  and  movements  were  analyzed  using 
conventional  discharge  histograms  (Fig.  5B)  and  raster  displays  (Fig.  5C).  In  addition,  ISI  rasters 
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were  plotted  that  displayed  the  time  of  occurrence  of  each  spike  on  the  x-axis  and  the  succeeding 
ISI  on  the  y-axis  (Fig.  5A).  To  analyze  vibration-related  changes  in  activity,  tte  occurrences  of 
individual  discharges  were  expressed  as  times  with  respect  to  vibration  onset  (Fig.  5  A  -  E,  left 
panels).  For  the  analysis  of  movement-related  changes,  spike  occurrences  were  expressed  wi 

respect  to  movement  onset  (Fig.  5 A  -  E,  right  panels).  .  .r/-.  j  *  -  a 

The  onsets  of  changes  in  neuronal  firing  rate  and  of  changes  in  EMG  were  dete^ned 
using  the  cumulative  sum  methods  (CUSUM)  (Ellaway,  1977;  Jiang  e  al  1991).  The  CUSUM 
at  a  given  time  is  the  total  number  of  discharges  accumulated  for  all  trials  from  some  starting  tmie. 
The  average  CUSUM  was  calculated  by  rescaling  this  count  by  dividing  it  by  the  number  of  trials 
(Fig.  5D).^The  CUSUM  for  an  impulse  train  for  which  the  probability  of  disdiarge  is  constant 
over  time  is  a  linearly  rising  curve  with  a  slope  equal  to  the  mean  firing  rate.  To  calculate  the  onse 
of  deviations  from  the  stationary  level  of  activity,  the  largest  epoch  of  linear  rise  in  the  aver  g 
CUSUM  prior  to  a  behaviorally  significant  event  was  labeled  by  visual  inspection.  A  linear  lea 
squares  interpolation  curve  was  calculated  for  this  period.  This  . 

tliough  the  epoch  containing  the  event.  The  standard  deviation  of  the  CUSUM  from  the  ihted 
curve  was  calculated  for  the  epoch  that  had  a  linear  rise.  The  computer  program  searched  forward 
in  time  to  find  the  first  change  in  the  CUSUM  from  the  curve  of  more  than  three  standard 
deviations  for  at  least  40  ms.  This  time  was  designated  as  the  onset  of  a  significant  change  in 

activity  jjjg  correlation  between  the  periodic  vibratory  stimulus  and  neuronal  discharges  was 
analyzed  using  a  phase  representation  method  (for  details,  see  Lebedev  et  al.,  1994).  The  phase  of 
each  spike  with  respect  to  the  stimulus  cycle  was  calculated  and  plotted  in  a  scattergram  as  a 
function  of  time  (see,  e.g..  Fig.  8B).  In  addition,  the  phase  relative  to  the  extrapolated  stimulus 
sinusoid  was  automatically  calculated  by  the  computer  for  the  epoch  preceding  vibration  onset. 

For  this  epoch,  phase  was  randomly  scattered.  Following  vibration  onset,  if  the  activity  was 
entrained  to  the  stimulus,  a  band  was  present  in  the  phase  scattergram  that  conesponded  to  the 
phase  of  preferential  response.  Cycle  distribution  histograms  also  were  calculated  that  represented 
the  probability  of  discharge  occurrence  in  relationship  to  the  phase  of  the  stimulus  cycle  (see,  e.g.. 
Fig.  8G). 

Statistical  Analyses 

The  characteristics  of  neuronal  activity  for  several  groups  of  neurons  (hping  specific  RF  types, 
located  in  given  cortical  areas,  and  so  on)  were  statistically  compared  using  a  multifactorial 
ANOVA  (with  the  Scheffe  post  hoc  test).  The  parametric  t-test  and  the  more  robust  nonparametric 
Mann-Whitney  U-test  were  used  for  two  group  comparisons. 


Results 

Cortical  Locations  and  Receptive  Fields 

Of  the  total  recorded  1,222  SI  neurons,  706  (-58%)  neurons  had  ^8  J^tes  of  more  than  10 
spikes/s  during  the  hold  period  of  the  paradigm.  Of  this  sample,  70/706  (~10  %)  neurons 
exhibited  bursty  firing  patterns  and  were  excluded.  Eight  of  the  remaining  636  neurons, 
were  excluded  from  consideration  because  they  had  either  a  small  number  of  ISIs  shorter  than  1ms 
(approximately  twice  the  spike  duration)  or  inconsistent  activity  patterns  during  the  experiment^ 
session.  Either  occurrence  could  indicate  inclusion  in  the  records  of  spikes  from  another  ne^by 
neuron.  Thus,  129/1222  neurons  (-11%  )  having  rhythmic  activity  patterns  were  selected  (Fig. 
3B)  The  distribution  of  surface  locations  of  the  recording  sites  and  their  location  in  representative 
sagittal  sections  are  shown  in  Fig.  2C.  No  clear  distribution  patterns  were  evident,  either  within  or 
across  cortical  areas  (areas  3a,  3b,  1,  and  2).  Ninety-one  neurons  were  tested  for  RFs.  This 
sample  contained  more  neurons  with  deep  RFs  than  those  with  cutaneous  RFs  or  no  cle^I^  (Fig. 
2B)  Cutaneous  RFs  most  frequently  were  associated  with  palmar  surface  of  the  hand.  Deep  Kr 
were  associated  with  movements  of  fingers  and  the  wrist  or  with  deep  hand  tissues. 
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Characteristics  of  Rhythmic  Activity 

The  activity  of  rhythmically  firing  neurons  that  occurred  during  the  hold  period  was  analyzed  to 
determine  its  frequency  characteristics  as  well  as  any  serial  dependencies  of  the  ISIs.  An  example 
of  spike  train  analyses  for  an  area  1  neuron  is  presented  in  Fig.  4.  The  activity  of  the  same  neuron 
during  task  execution  is  illustrated  in  Fig.  5.  In  Fig.  4A,  a  scattergram  of  ISIs  during  a  500-ms 
epoch  preceding  vibration  onset  is  plotted  for  40  consecutive  trials.  This  scattergram  shows  that 
the  ISI  distribution  remained  virtually  unchanged  from  trial  to  trial.  In  Fig.  4B,  a  joint  interval 
scattergram  is  presented.  Fig.  4C,  D,  and  E  display  ISI  distribution,  ED,  and  RD  histograms, 
respectively.  During  the  hold  period,  this  neuron  was  rhythmically  active  at  ~39  Hz.  The  majority 
of  ISIs  were  distributed  within  the  range  of  +  5  ms  around  the  rhythmic  period  (-25.6  ms).  Note, 
however,  that  a  small  number  of  outlying  short  and  long  ISIs  occurred.  ED  and  RD  histograms 
were  not  substantially  different. 

Statistical  analyses  of  the  first  normalized  ED  peaks  equation  1)  and  the  rhythmic 
frequencies  (/ ;  equation  2)  for  the  sample  of  rhythmically  active  neurons  did  not  show  any 
significant  differences  depending  on  cortical  location  of  the  neurons  nor  on  their  RF  type.  The 
means  and  the  standard  deviations  were,  for  £"/,  2.62  ±  0.78  (Fig.  3B)  and,  for  /,  32.1  +  5.4 

Hz  (Fig.  3E).  .  ,  r  u  u 

We  observed  two  features  of  spike  trains  that  seem  unlikely  for  a  renewal  model  of  rhythm 

generation  (Perkel  et  al.,  1967).  These  features  were  (1)  multimodal  ISI  distributions  and  (2)  the 
occurrence  of  a  small  but  noticeable  number  of  ISIs  at  less  than  the  modal  interval  that  formed 
diagonal  bands  in  joint  interval  scattergrams.  Multimodal  ISI  distributions  were  obsen^ed  for 
31/129  neurons  (-24%).  An  example  of  an  area  1  neuron  with  a  multimodal  ISI  distribution  is 
presented  in  Fig.  6A-D  (also  see  Figs.  9  and  1 1).  The  ISI  distribution  of  this  neuron  contained 
peaks  at  the  modal  interval  (7)  and  at  twice  that  interval  (27,  Fig.  6B).  The  joint  interval 
scattergram  for  this  neuron  had  four  clusters  of  ISIs  pairs  around  points  (7,  7),  (7,  27),  (27, 

T),  and  (27, 27).  ED  and  RD  histograms  were  approximately  the  same  for  this  neuron.  For 
32/129  neurons  (-25%),  we  observed  the  occurrence  of  ISIs  that  were  shorter  than  modal  interval. 
Six  neurons  (-5%)  exhibited  both  these  short  ISIs  and  multimodal  ISIs.  Short  ISIs  that  form 
diagonal  bands  in  joint  interval  plots  have  been  attributed  to  “interrupting  spikes”  (Surmeier  and 
Towe,  1987a,  1987b;  Siebler  et  al.,  1991).  Joint  interval  scattergrams  for  such  spike  trains 
contain  diagonal  bands  that  connect  points  (0,  7)  and  (7,  0).  These  bands  occur  when  the  sum 
of  the  ISIs  that  immediately  precede  and  succeed  the  intermpting  spikes  is  equal  to  the  modal 
interval: 

Tj  +  T2  =  T,  (3) 

where  7;  is  the  ISI  preceding  an  interrupting  spike  and  T2  is  the  succeeding  ISI.  Records  for  an 
area  3b  neuron  that  exhibited  interrupting  spikes  are  presented  in  Fig.  6E  -  H.  In  addition  to  a 
cluster  of  ISIs  around  7,  shorter  ISIs  were  present  (from  2  ms  to  7,  Fig.  6F).  In  the  joint 
interval  scattergram,  the  presence  of  short  intervals  resulted  in  a  diagonal  band.  Note  that  peaks  in 
the  RD  histogram  were  less  pronounced  than  those  in  the  ED  histogram. 

The  comparison  of  ED  and  RD  histograms  for  the  total  sample  of  rhythmically  firing  cells 
showed  a  tendency  for  peaks  in  RD  histograms  to  be  smaller  in  amplitude  than  those  in  ED 
histograms.  The  results  of  a  factorial  ANOVA  indicated  that,  for  the  neurons  having  interrupting 
spikes,  the  differences  between  ED  and  RD  peaks  were  significantly  greater  compared  with  the  rest 
of  the  total  sample.  For  these  neurons,  the  means  and  standard  deviations  of  the  differences 
between  the  first  two  ED  and  RD  peaks  were  13.4  +  9.7%  and  13.8  +  9.2%,  respectively  (p  < 
0.0002;  paired  r-test).  When  the  neurons  with  interrupting  spikes  were  excluded  from  the  total 
sample,  a  small  but  statistically  significant  difference  was  found  only  for  the  first  peaks  (1.8  + 
5.0%;  p  <  0.0002). 

Vibration-Related  Activity 

Changes  in  the  activity  of  SI  neurons  in  this  experimental  paradigm  commonly  occurred  at 
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vibratory  go-cue  onset  and  often  preceded  movement  onset  (Nelson,  1988).  We  an^yzed  the 
changes  in  activity  of  rhythmically  firing  neurons  that  occurred  in  association  with  these 
behaviorally  significant  events.  For  67/129  rhythimcally  active  neurons  (-52%),  the  firing  rate 
changed  with  vibration  onset.  Decreases  in  activity  were  the  most  common  responses  to  vibration 
(6Qlil  cases  -89%).  An  example  of  a  typical  response  to  vibration  exhibited  by  an  area  1  neuron 
is  oresented  in  Fig.  7 A.  The  neuron’s  firing  rate  decreased  following  palinar  vibration  with  a 
latency  of  -28  ms  The  vibration-related  activity  decrease  consisted  of  both  transient  and  sustained 
changL  in  activity  (schematic  illustration.  Fig.  7E).  The  analyses  of  the  latency  of  vibratory 
respo^nses  are  presented  in  Fig.  7F.  This  latency  (-26  ms  at  57  Hz)  was  not  significantly  different 

as  a  function  of  cortical  location  or  RF  type. 

In  some  cases,  neuronal  activity  during  ongoing  vibration,  although  decreas^,  was 
stimulus-entrained.  An  example  of  an  area  3a  neuron  with  activity  entrained  to  a  27  Hz  vibratory 
stimulus  is  presented  in  Fig.  8.  During  the  hold  period,  the  neuron  was  rhythmically  active  at  -30 
Hz.  A  transient  decrease  in  activity  occurred  at  -27  ms  after  vibration  onset  and  lasted  for  -50  ms 
IFig.  8A).  Then,  the  activity  recovered,  although  it  was  at  a  lower  level  thM  that  during  the  ho 
period.  Moreover,  neuronal  discharges  were  entrained  to  the  ongoing  27  Hz  vibration.  The 
entrainment  is  clearly  seen  in  phase  raster  (Fig.  8B)  and  in  the  cycle  distribution  histogram  (  ig. 

8G)  The  ISI  distribution  for  the  epoch  of  vibratory  stimulation  contained  peaks  at  the  vibrato^ 
period  and  at  twice  that  period  (Fig.  8F).  Generally,  entrainment  was  obse^ed  more  often  when 
stimulating  at  27  Hz,  which  was  close  to  the  average  rhythrmc  frequency  of  the  studied  neurons. 

Of  4 1  neurons  recorded  during  stimulating  at  this  frequency,  12  neurons  (-29%)  exhibited 
vibration-entrained  activity.  The  majority  of  the  data  were  collected  while  stimulating  at  57  Hz.  At 
this  frequency,  12/129  neurons  (-9%)  exhibited  vibratory  entrainment.  At  the  highest  stimulus 
frequency  (127  Hz),  vibratory  entrainment  rarely  was  observed  (1/38;  -3%). 

Premovement  Activity 

The  activity  of  each  rhythmically  firing  neuron  changed  prior  to  movement  onset  for  at  least  one 
movement  direction.  The  earliest  change  of  firing  rate  from  the  stabilized 

vibratory  stimulation  was  designated  as  the  onset  of  premovement  activity  (PMA)  (Nelson,  1988). 
PMA  was  detected  in  124/129  instances  for  flexion  movements  and  m  1 1^/129  instances  for 
extension  movements.  Decrease  in  firing  rate  was  the  most  frequent  type  of  PMA  (-76%  of  PMA 
cases-  Fig.  9F).  Figure  9A  illustrates  this  pattern.  Records  for  an  area  2  neuron  are  presented. 

The  firing  rate  of  this  neuron  did  not  change  after  vibration  onset.  However,  its  activity  was 
dramatically  decreased  at  -40  ms  prior  to  the  onset  of  extension  movements.  The  PMA  onsets 
with  respect  to  movement  onset  were  not  statistically  different  as  a  function  of  corticm  location  or 
RF  type  These  onsets  also  did  not  differ  between  instances  of  activity  increases  and  decreases. 

To  estimate  the  temporal  relationship  between  PMA  and  the  earliest  movement-related 
peripheral  afferent  signals,  PMA  onsets  were  compared  with  EMG  onsets.  ^he  EMG  onsets  were 
analyzed  for  several  forearm  and  arm  muscles  (Lebedev  et  al.,  1994;  Nelson,  1987,  Nelson  et  al., 
1991).  A  minimal  afferent  peripheral  conduction  time  of  1 1  ms  was  subtracted 
EMG  onsets  (Wiesendanger  and  Miles,  1982)  to  yield  an  estimate  of  the  tiine  at  which  afferent 
signals  resulting  from  muscle  activity  might  reach  the  cortex.  The  earliest  tm^  at  which  this  might 
oa:ur  was  thus  100  ms  prior  to  movement  onset,  whereas  given  the  average  EMG  onset,  it  was  60 
ms  prior  to  movement  onset.  A  substantial  number  of  PMA 

occurred  after  the  earliest  estimated  time  of  afferent  input  arrival  (194/242,  ^ 

PMA  onsets  occurred  after  the  average  estimated  time  of  afferent  input  arrival  (133/242,  f’)- 

During  activity  decreases,  rhythmic  activity  often  was  dismpted  rather  than  modulatea 
(e  g  Fig  9A).  In  the  cases  when  the  rhythmic  pattern  of  discharges  was  presented,  the  rhythme 
frequency  decreased  slightly  (by  -20%;  e.g.,  sustained  vibratory  response  m  Fig.  5A).  Activity 
increases  usually  occurred  during  movements  and  were  charactenzed  by  transitions  from  rhythrmc 
to  nonrhythmic  firing.  In  Fig.  10,  records  for  an  area  1  neuron  are  presented  that  show  activation 
during  voluntary  extensions.  This  neuron  also  was  activated  by  passive  extension  of  the  second 
digit  at  the  metacarpophalangeal  joint  (Fig.  lOB).  Rhythmic  firing  decreased  at  -27  ms  following 
vibration  onset  and,  additionally,  at  -34  ms  prior  to  extension  movement  onset.  Then,  a 
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pronounced  activation  followed  at  ~45  ms  after  movement  onset.  This  activation  was  characterized 
by  a  qualitative  change  in  the  ISI  distribution  (Fig.  lOD  and  E).  The  ISI  distribution  shifted  to  the 
left,  indicative  of  the  faster  firing  rate.  However,  several  ISIs  occupd  at  the  minimum  value  of 
1.2'ms,  possibly  indicating  the  presence  of  small  bursts  of  activity  in  this  otherwise  rhythmically 
spiking  neuron.  In  some  instances,  transitions  to  irregular  ISI  patterns  occurred  without 
substantial  changes  in  the  mean  firing  rate  (e.g.,  -200  ms  after  movement  onset  in  Fig.  6E  and 
1  lA). 

Activity  patterns  were  analyzed  to  determine  if  PMA  types  were  dependent  upon  the 
direction  of  subsequent  movement.  If  the  sign  of  the  activity  change  was  opposite  for  flexion  and 
extension  or  if  PMA  occurred  only  for  one  movement  direction,  this  instance  of  PMA  was 
classified  as  directional.  If  the  activity  change  was  of  the  same  sign  for  both  movement  directions, 
the  PMA  was  classified  as  nondirectional.  Nondirectional  PMA  occurred  most  frequently  (84/129 
neurons,  -65%);  directional  PMA  was  observed  for  45/129  neurons  (-35%).  An  example  of  an 
area  3a  neuron  with  a  nondirectional  PMA  pattern  is  presented  in  Fig.  11.  For  both  flexions  (Fig. 

1  lA)  and  extensions  (Fig.  1  ID),  premovement  decreases  in  activity  of  this  neuron  occurred  at  -70 
ms  prior  to  movement  onset. 

Comparison  with  Other  SI  Neurons 

We  compared  the  activity  pattern  exhibited  by  rhythmically  firing  neurons  with  the  types  of  activity 
of  SI  neurons  that  have  been  documented  previously  for  this  experimental  paradigm  (Nelson, 

1988;  Nelson  et  al.,  1991).  The  typical  activity  pattern  of  rhythmically  firing  neurons  resembled  a 
mirror  image  of  the  activity  of  quickly  adapting  (QA)  neurons  recorded  during  similar  behaviors 
(Nelson  et  al.,  1991).  QA  neurons  were  activated  at  about  the  same  time  following  vibratory  cue 
onset  (-30  ms)  and  prior  to  movement  (-80  ms)  as  the  firing  rate  of  rhythmically  firing  SI  neurons 
decreased. 

Discussion 

Activity  of  Rhythmically  Firing  Neurons  During  Trained  Motor  Tasks 

We  examined  the  activity  of  a  population  of  rhythmically  firing  neurons  in  monkey  SI  while  the 
animals  performed  trained  motor  tasks.  The  changes  in  activity  of  these  neurons  had  common 
features  during  the  initiation  of  vibratory-cued  movements.  Approximately  one-half  of 
rhythmically  firing  neurons  responded  to  vibratory  go-cues,  typically  with  decreases  in  firing  rates 
at  -25  ms  after  vibration  onset.  In  some  cases,  although  the  firing  rate  was  decreased  somewhat, 
the  activity  was  entrained  to  the  frequency  of  the  ongoing  vibration.  It  is  possible  that  the 
entrainment  observed  was  caused  by  IPSPs  rather  than  by  EPSPs  (Cowan  and  Wilson,  1994; 
Lytton  and  Sejnowski,  1991;  Stern  and  Wilson,  1994).  The  best  vibratory  frequency  for  this 
entrainment  was  close  to  the  population's  average  frequency  of  rhythmic  activity  (-30  Hz).  This 
observation  is  consistent  with  a  model  of  texture  perception  proposed  by  Ahissar  and  Vaadia 
(1990).  In  this  model,  cortical  oscillators  are  key  elements  in  neuronal  circuits  that  analyze  the 
temporal  properties  of  somatosensory  discharges  using  a  phase-locked  loop  algorithm  (Horowitz 
and  Hill,  1980).  The  model  predicts  that  oscillatory  neurons  would  follow  somatosensory  stimuli 
if  the  stimulus  frequency  is  close  to  the  frequency  of  spontaneous  oscillations  of  these  neurons. 
Although  our  experiments  were  not  specifically  designed  to  examine  the  role  of  rhythmically  active 
neurons  in  somatosensory  perception,  our  results  generally  support  this  prediction. 

Decreases  in  the  activity  of  rhythmically  firing  neurons  often  occurred  at  -57  ms  prior  to 
movement  onset.  Thus,  a  typical  activity  pattern  consisted  of  a  decrease  in  firing  rate  after 
vibration  onset  and  an  additional  decrease  in  firing  rate  prior  to  movement  onset.  This  activity 
pattern  is  complementary  to  the  pattern  previously  reported  for  QA  neurons  that  were  recorded 
during  the  performance  of  this  behavioral  task  (Nelson  et  al.,  1991).  QA  neurons  exhibit  a  short- 
duration  burst  of  spikes  after  vibration  onset.  Then  their  activity  decreases  and  is  reactivated  prior 
to  movement  onset.  Nelson  et  al.  (1991)  suggested  that  QA  activity  patterns  may  occur  because  of 
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gating  of  SI  sensory  responsiveness  at  behaviorally  significant  times.  The  function  of  gating  may 
be  to  enhance  sensory  inputs  that  are  important  for  the  current  behavioral  task  while  suppressing 
Xs  am  not  (Ch^n  and  Woodward,  1982;  Chapman  et  al.  1988;  Coquery,  1978;  Dyhre- 
Poulsen  1978;  Nelson,  1988;  Rushton  et  al.,  1981).  Because  of  the  correspondence  between  the 
activity  patterns  of  rhythmically  firing  and  QA  neurons,  it  seems  reasonable  to  assume  that 
rhythmically  firing  neurons  may  participate  in  gating  SI  activity. 

^  Centrally  generated  as  well  as  peripheral  inputs  to  SI  neurons  may  be  involved  m 
premovement  gating  of  somatosensory  activity  (Chapin  and  Woodw^d,  1982;  Ch^man  et  a  ., 
1988‘  Coquery,  1978;  Dyhre-Poulsen,  1978;  Lebedev  et  al.,  1994;  h^lson,  1987;  Nelson  et  a  ., 
1991-  Soso  and  Fetz,  1980).  Changes  in  the  activity  of  rhythnucally  firing  neurons  that  occurred 
before  movement  usually  happened  after  the  earliest  EMG  onset.  Therefore,  these  chanps 
certainly  could  be  related  to  peripheral  afferent  signals  associated  with  the  onset  of  muscle  activity 
that  preceded  movement  onset  as  detected  by  the  change  m  handle  position.  However,  most 
premovement  activity  changes  were  of  the  same  sign  for  both  fle^on  and  extension  trials  (*at  , 
Siey  were  nondirectional).  Peripheral  reafferent  signals  commonly 

of  proprioceptive  origin  (Cohen  et  al.,  1994;  Soso  and  Fetz,  1980).  Many  of  the  rhythnucally 
firing  SI  neurons  in  our  sample  had  deep  RFs.  Nonetheless,  they  usually  exlubited  nondirectiona 
nremovement  activity  patterns.  Therefore,  the  premovement  changes  m  activity  of  rhythmcally 
firing  SI  neurons  probably  are  not  simply  replicas  of  peripheral  reafferent  signals  but  rather  mig 
reflect  intracortical  processing  of  peripheral  information. 


Rhythmically  Active  SI  Neurons  and  LFP  Oscillations 

Our  present  observations  suggest  that  there  are  similar  features  between  the  activity  patterns  of 
rhythmically  firing  SI  neurons  and  the  LFP  oscillations  in  sensonmotor  cortex  reported  by  others. 
The  mean  firing  frequency  of  the  observed  rhythmically  firing  SI  neurons  (~30  Hz)  was  close  to 
the  frequency  often  described  for  LFP  oscillations  (Bouyer  et  al.,  1981,  1987;  Murthy  and  Fetz 
1992'  Rougeul  et  al.,  1979;  Sanes  and  Donoghue,  1993).  In  a  manner  similar  to  LFP  oscillations, 
rhvthmic  activity  of  single  SI  neurons  was  disrupted  with  somatosensory  stimulation  and/or 
voluntary  movements.  LFP  oscillations  occur  in  sensorimotor  cortex  in  the  cat  (Bouyer  etaL 
198 1  1987)  and  the  monkey  (Rougeul  et  al.,  1979)  during  focused  attention  and  immobility  but 
usually  disappear  during  movements.  Disruptions  of  LFP  oscillations  m  monkey  motor  cortex 
also  occur  prior  to  movements  performed  during  trained  motor  tasks  (Sanes  and  Donoghue,  1993). 
In  addition,  Murthy  and  Fetz  (1992)  observed  that  LFP  oscillations  m  monkey  motor  and 
somatosensory  cortices  occurred  during  exploratory  movements  requiring  directed  attention. 

Thus  LFP  oscillations  may  vary  as  a  function  of  attentive  behavior  and  task  retirements  (tor 
review  see  Gray  1994).  Experimental  conditions  of  Sanes  and  Donoghue  (1993)  are  the  most 
similar’to  ours.  Disruptions  of  LFP  oscillations  in  their  experiments  occurred  ^ound  the  same 
time  as  the  disruptions  of  rhythmic  firing  of  individual  SI  neurons  in  our  experiments.  We 
suggest,  therefore,  that  LFP  oscillations  and  rhythmic  firing  of  SI  neurons  may  be  related. 

Some  characteristics  of  rhythmic  activity  of  individual  neurons,  however,  differ  from  those 
reported  for  LFP  oscillations.  LFP  oscillations  usually  occur  as  occasional  rhythimc  bursts,  lUU  - 
200  ms  in  duration  and  often  with  a  variable  frequency  from  one  burst  to  another  (Ec^om  et  al., 
1988'  Gray  and  Singer,  1989;  Murthy  and  Fetz,  1992;  Sanes  and  Donoghue,  1993,  Singer,  1993). 
The  activity  patterns  of  individual  rhythmically  firing  SI  neurons  were  consistent  from  trial  to  trial 
and  were  only  briefly  interrupted  during  vibratory  stimulation  and  movements.  We  suggest  that 
individual  rhythmically  firing  SI  neurons  may  maintain  relatively  independent  activity  Further,  u 
is  possible  that  these  neurons  may  intermittently  evoke  oscillations  in  other  neurons.  Murthy  an 
Fetz  (1994)  observed  sensorimotor  cortical  neurons  that  intermittently  exhibited  episodes  or 
rhvthmic  activity  during  episodes  of  LFP  oscillations.  These  neurons  could  be  driven  by  neurons 
with  sustained  rhythmic  activity.  Driven  oscillations  of  cortical  neurons  inay 
through  cortico-cortical  and  thalamocortical  interactions  (Limas,  1990;  Limas  et  al.,  1991).  thus, 
rhythmically  firing  neurons  may  act  to  regulate  cortical  population  oscillations. 
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Possible  Models  of  Rhythmic  Activity 


Mechanisms  causing  rhythmic  activity  of  SI  neurons  remain  unclear.  One  possibility  is  that 
rhythmically  active  cells  are  integrate-and-fire  neurons,  and  their  regular  firing  occurs  because  of 
stable  synaptic  input  (Segundo  et  ah,  1968).  Regular  spiking  cortical  pyramidal  neurons  discharge 
very  rhythmically  in  an  integrate-and-fire  mode,  when  they  receive  steady  inputs  (McCormick  et 
al.  1985-  Softky  and  Koch,  1993).  Integrate-and-fire  neurons  respond  to  changes  in  the  intensity 
of  synaptic  input  by  changing  their  firing  frequency.  If  rhythmically  firing  SI  neurons  are 
integrate-and-fire  neurons  activated  by  stable  peripheral  inputs,  then  their  firing  frequency  should 
follow  changes  in  the  intensity  of  peripheral  inputs,  for  example,  as  might  occur  cmring  ,  ,  , 

movements.  However,  the  rhythmic  firing  of  SI  neurons  is  disrupted  rather  than  t^ing  modulated 
during  movements.  Thus,  modeling  rhythmically  firing  neurons  as  integrate-and-fire  neurons  does 

not  describe  all  of  our  observations.  .u  •  j  •  • 

Alternatively,  rhythmic  firing  may  result  from  extrinsic  oscillatory  drive.  Rhythmic  driving 

of  cortical  neurons  produced  by  intracortical  and/or  subcortical  inputs  previously  has  been 
proposed.  Chose  and  Freeman  (1992)  suggested  that  visual  cortical  neurons  may  be  driven  by 
rhythmically  active  cells  located  in  the  lateral  geniculate  nucleus.  Jagadeesh  et  al.  (1992) 
demonstrated  rhythmic  excitatory  inputs  to  visual  cortical  neurons  using  in  vivo  patch  clamp 
recording  techniques.  The  source  of  rhythmic  drive  may  be  from  inhibitory,  rather  than  excitatory, 
neurons  (Cowan  and  Wilson,  1994;  Llinas,  1990;  Llinas  et  al.,  1991).  Some  of  our  observations 
are  consistent  with  an  extrinsic  drive  model.  We  observed  ISIs  for  some  neurons  that  occurred  at 
multiples  of  the  modal  interval.  This  may  have  occurred  because  of  a  neuron’s  occasional  failures 
to  respond  to  driving  inputs  as  previously  suggested  (Ahissar  and  Vaadia,  1990).  In  addition,  we 
observed  interrupting  spikes  that  may  reflect  spikes  that  are  uncorrelated  with  the  rhythmic  drive 
(Surmeier  and  Towe,  1987a,  1987b).  However,  from  the  manner  in  which  these  neurons  were 
recorded  we  cannot  absolutely  eliminate  the  possibility  that  interrupting  spikes  or  ISIs  at  multiples 
of  the  modal  interval  could  reflect  occasional  spurious  inclusion  spikes  from  adjacent  neurons  with 
similar  waveforms  or  exclusion  of  spikes  from  the  neuron  in  question.  Further  experiments  are 


required  to  clarify  this  issue.  •  .  i 

Rhythmic  driving  of  SI  neurons  could  arise  from  regular  discharges  of  peripheral  atterents. 
Single  peripheral  afferent  fibers  may  have  very  regular  firing  patterns  (e.g.,  Burke  et  al,  1987). 
However  the  frequency  of  rhythmic  discharge  of  peripheral  afferents  is  modulated  during 
movements  (Burke  et  al,  1987;  Prochazka,  1985).  Since  rhythmic  activity  of  SI  neurons  typically 
was  disrupted  periodically,  these  neurons  were  unlikely  to  be  driven  by  peripheral  afferent 
discharges  throughout  all  phases  of  our  behavioral  task.  Moreover,  the  SI  neurons  that  seem  to 
have  been  driven  most  securely  by  peripheral  inputs  did  not  exhibit  rhythmic  activity.  These 
neurons  which  had  cutaneous  or  deep  RFs  faithfully  responded  to  vibratory  go-cues  with  entrained 
discharges  (Lebedev  et  al.,  1994).  However,  they  did  not  fire  rhythmically  during  the  hold  period 
of  the  paradigm  (unpublished  observations).  Therefore,  it  seems  unlikely  that  peripheral  drive  is 
the  major  source  of  rhythmic  firing  of  SI  neurons  under  these  experimental  conditions. 

Rhythmic  activity  of  some  of  the  studied  neurons  could  be  the  result  of  intrinsic  oscillatory 
properties  rather  than  extrinsic  factors.  For  many  rhythmically  firing  neurons  (72/129;  ~56%)  we 
did  not  observe  signs  of  extrinsic  drive  (that  is,  interrupting  spikes  and  multimodal  distributions  ot 
the  ISIs)  Moreover,  ED  and  RD  histograms  for  these  neurons  were  virtually  identical. 
Mountcastle  et  al.  (1969, 1990)  reported  that  ED  and  RD  histograms  were  dramatically  different 
when  the  firing  of  SI  neurons  was  driven  by  peripheral  vibration.  Absence  of  such  differences 
between  ED  and  RD  histograms  for  many  rhythmically  active  SI  neurons  may  indicate  either  that 
these  neurons  were  not  extrinsically  driven  or  that  they  were  driven  by  very  secure  oscillatory 


inputs.  _  - 

Intrinsically  oscillating  cortical  neurons  have  been  reported.  Two  types  oi  neurons 
recorded  in  vitro  have  been  demonstrated  to  intrinsically  oscillate  at  10  -  50  Hz  within  4  oi 
guinea  pig  frontal  cortex  (Llinas  et  ah,  1991).  Neurons  of  the  first  type  -  that  is,  broadTrequency 
oscillators  -  increased  their  firing  frequency  with  membrane  depolarization.  Neurons  of  the  second 
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type  the  narrow-frequency  oscillators  showed  little  change  in  firing  frequency  as  a  function  of  the 
level  of  depolarization.  The  narrow-frequency  oscillatory  neurons  have  been  identified 
morphologically  as  sparsely  spinous  interneurons  having  axon  collaterals  in  layers  3  and  4  and  are 
probably  inhibitory  intemeurons.  Pyramidal  cortical  neurons  typically  do  not  generate  fast  (20  -  50 
Hz)  oscillations  intrinsically  (SiWa  et  al.,  1991).  Rather,  their  actiyity  is  thought  to  be  moddat^ 
by  rhythmic  IPSPs  probably  generated  by  cortical  interneurons  (Cowan  and  Wilson,  1994,  Stem 
and  Wilson,  1994;  but  see  also  Nunez  et  al.,1992).  Thus  observations  tend  to  indicate  that 
rhythmically  firing  neurons  in  the  cortex  may  be  inhibitory  intemeurons. 


Rhythmically  Active  Neurons  and  Tonic  Inhibition  of  SI  Activity 


Since  cortical  oscillatory  neurons  are  likely  to  be  inhibitory  intemeurons  (Cowan  and  Wilson, 
1994-  Llinas  et  al.,  1991;  Stern  and  Wilson,  1994),  there  is  the  intriguing  possibiluy  that 
rhvthmically  active  SI  neurons  may  be  involved  in  tonic  inhibition  of  SI  actiyity.  The  role  or 
inhibitory  intemeurons  as  essential  regulators  of  cortical  activity  is  suggested  by  the  obse^ation 
that  -70%  of  intracortical  connections  are  inhibitory  (White,  1989).  Tonic  inhibition  produced  by 
local  intemeurons  has  been  suggested  to  control  RF  size  of  cortical  neurons  and  to  increase  spatial 
and  temporal  contrast  (Brooks,  1959;  Dykes  et  al.,  1984).  We  further  suggest  that  inhibUory 
intemeurons  may  gate  SI  activity  during  behavior.  The  activity  patterns  of  rhythmically  finng 
neurons  are  complimentary  to  the  activity  patterns  of  quickly  adapting  (QA)  neurons  recorded 
during  this  behavior  (Nelson  et  al.,  1991).  This  observation  is  similar  to  that  of  Wilson  et  al. 
(1994)  who  observed  that  putative  GABA-ergic  intemeurons  and  pyramidal  neurons  exhibited 
complementary  patterns  of  activity  while  the  monkeys  performed  visual  and  oculomotor  tasks. 
These  authors  suggested  that  intemeurons  may  hyperpolarize  pyramidal  neurons,  whereas 
pyramidal  neurons  may  depolarize  intemeurons.  Similar  reciprocal  connections  may  be  present 

between  rhythmically  firing  and  QA  SI  neurons.  ,  , .  •  i 

Rhythmically  firing  neurons  are  candidates  for  one  element  of  a  model  that  was  previously 

proposed  to  account  for  observed  patterns  of  activity  of  QA  SI  neurons  (Nelson  et  ^.,  1991). 
task  phase  element  was  proposed  that  may  suppress  peripheral  and  motor  inputs  to  QA  neurons. 
In  Fig.  12,  a  model  is  presented  that  is  a  modified  form  of  the  previous  one.  This  model  includes 
output  neurons  (such  as  corticocortical,  corticothalamic  or  corticofugal),  input  neurons 
(presumably,  layer  IV  spiny  cells)  (Jones,  1975)  and  two  types  of  inhibitory  intemeurons, 
rhythmic  and  phasic  intemeurons.  This  model  assumes  the  following  functional  relationships. 
Rhythmic  intemeurons  exhibit  sustained  firing  patterns,  whereas  phasic  inhibitory  interneurons 
respond  to  inputs  with  transient  bursts  of  spikes.  The  model  includes  neurons  that  receive  direct 
excitatory  inputs  from  other  regions.  The  model  has  been  conceived  such  that  rhythmically  tiring 
neurons  tonically  inhibit  output  neurons.  Output  neurons  may  be  disinhibited  becmse  their 
inhibitors  themselves  receive  inhibitory  inputs  from  phasic  intemeurons.  During  increases  in  ol 
input,  phasic  inhibitory  neurons  would  be  transiently  activated.  Consequently,  by  this  model, 
rhythmically  firing  neurons  would  be  transiently  inhibited.  This  scheme  suggests^  that  output 
neurons  generate  bursts  of  spikes  when  released  from  tonic  inhibition  and  are  excited  by  input 
signals.  Quickly  adapting  activity  patterns  of  output  neurons  may  be  important  for  signalling 


transient  sensory  events.  . 

The  model  shown  in  Fig.  12  is  consistent  with  the  data  on  intracortical  connectivity 
(Chagnac-Amitai  and  Connors,  1989;  Istvan  and  Zarzecld,  1994;  Jones,  1975).  Nonpyramidd 
inhibitory  neurons  that  may  regulate  activity  of  other  cortical  neurons  previously  have  been  stuched 
(Kawaguchi,  1993;  Kawaguchi  and  Kubota,  1993;  Jones,  1975).  One  type  of  neuron,  the  double 
bouquet  cell  (Ramon  y  Cajal,  191 1),  has  been  shown  to  project  to  nonpyramidal  cortical  cells 
(Somogyi  and  Cowey,  1981).  These  cells  are  characterized  by  vertically  oriented  axonal  arbors. 
Recent  findings  suggest  that  GABA-ergic  cortical  neurons  with  this  type  of  axonal  arbors  ^e 
immunoreactive  for  calbindin-28,  whereas  GABA-ergic  neurons  with  horizontally  oriented  axonal 
arbors  (for  example,  the  basket  cells)  are  immunoreactive  for  parvalbumin  (Hendry  et  al.,  1989; 
Kawaguchi  and  Kubota,  1993).  Moreover,  calbindin-28-iminunoreactive  neurons  respond  to 
current  injections  by  firing  phasically,  whereas  parvalbumin-immunoreactive  neurons  respond  by 
firing  repetitively  with  little  or  no  spike  frequency  adaptation  (Kawaguchi  and  Kubota,  1993). 
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Therefore,  c3lbinclin-28  cind  p&rvslbuirun-immunorecictive  neurons  sre  likely  cundidutes  for  the 
phasic  and  rhythmic  inhibitory  elements  of  the  model,  respectively. 

This  model's  circuitry  may  be  capable  of  initiating  network  oscillations  following  the 
scheme  proposed  by  Llinasetal.  (1991).  Rhythmically  firing  inhibitory  neurons  may  induce 
rhythmic  activity  in  output  neurons.  Output  neurons,  in  turn,  may  induce  rhythmic  activity  in  their 
target  regions,  for  example  in  the  thalamus  or  in  other  cortical  areas.  Regions  projecting  back  to 
the  cortex  would,  thus,  close  the  loop  which  would  be  necessary  for  the  maintenance  of  network 
oscillations. 
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Figure  Legends 

Fie  1  Schematics  of  the  manipulandum,  the  feedback  display,  and  the  task.  A:  Apparatus  for 
studying  wrist  movements.  B:  Visual  display  of  wrist  position  that  consisted  of  light-emitting 
diodes.  C:  Diagram  of  experimental  paradigm. 

Fie.  2.  Cortical  locations  of  rhythmically  firing  SI  neurons.  A:  A  drawing  of  the  dorsolateral 

view  of  the  brain.  Central  sulcus  (CS)  and  intraparietal  sulcus  (IPS)  are  indicated.  B:  A  table 

showing  the  number  of  rhythmically  firing  neurons  by  receptive  field  (I^)  type  recorded  m  each 
SI  region  (areas  3a,  3b,  1,  and  2).  C:  Locations  of  recording  sites  for  six  monk^s.  The 
illustration  for  each  monkey  consists  of  four  panels.  The  top  panel  presents  a  surface  map  o 
electrode  penetrations  with  respect  to  CS  and  IPS.  The  penetrations  where  rhythimcally  firing 
neurons  were  recorded  are  indicated  by  larger  marks.  The  lower  panels  illustrate  three  sagittal 
sections  through  the  cortex  (lateral,  intermediate,  and  medial).  Locations  of  the  neurons  are  shown 
as  projections  to  the  nearest  of  the  three  sections.  RF  types  are  marked  according  to  the  convention 

of  panel  B. 

Fig.  3.  Selection  of  the  group  of  rhythmically  firing  neurons.  A:  A  computer-generated 
of  an  expectation  density  (ED)  histogram.  The  spike  tram  constructed  to  have  a  no^al  distribution 
of  interspike  intervals  corresponding  to  30  Hz  and  a  coefficient  of  variation  of  0. 15.  R  is  th 
stabilized  level  of  the  histogram,  Ej,  E2  are  the  heights,  and  t;,  ^2  are  the  times  of  the  initia 
peaks.  B:  Histograms  of  the  distribution  of  the  ratio  of  the  first  peak  to  stabilized  ED  level  (Ly  / 

R)  for  rhythmically  firing  and  nonrhythmically  firing  neurons.  C-D:  scattergrams  illustratmg 
selection  of  rhythmically  firing  neurons  based  on  the  heights  of  ED  peaks  . 

interpeak  inteiwal.  Jitter  is  the  time  of  the  first  peak  minus  the  time  between  the  first  and  second 
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peaks  divided  by  the  time  of  the  first  peak.  E;  Histogram  of  the  distribution  of  the  firing  frequency 

for  the  total  sample  of  rhythmically  firing  SI  neurons. 

Fig.  4.  Analyses  of  rhythmic  activity  (~39  Hz;  same  records  as  in  Fig.  5).  Spike  trains  for  the 
hold  period  (500  ms  epoch  preceding  vibration  onset)  were  analyzed.  A:  Hold  period  interspike 
interval  (ISI)  scattergram  for  40  consecutive  trials.  Each  column  represents  a  trial.  B:  Joint 
interval  scattergram.  C:  ISI  distribution  histogram.  D:  Expectation  density  histogram.  E:  Renewal 
density  histogram.  In  panels  C  -  E,  bin  width  is  1  ms. 

Fig.  5.  Example  of  changes  in  neuronal  activity  associated  with  57  Hz  vibratory  go-cue 
presentation  and  movement.  Records  in  the  left  parts  of  panels  A  -  E  are  centered  on  vibration 
onset  whereas  those  in  the  right  parts  are  centered  on  movement  onset.  The  activity  pattern  of  a 
rhythmically  firing  (~39  Hz)  area  1  neuron  exhibited  during  flexion  trials  is  illustrated.  A:  Raster 
displays  of  interspike  intervals  (ISIs).  For  each  spike,  these  displays  show  the  time  of  its 
occurrence  on  the  x-axis  and  plot  the  next  ISI  on  the  y-axis.  B:  Histograms  of  discharge  rate  (bin 
width  =  5  ms).  Dotted  lines  represent  average  movement  onset  (left  panel)  and  average  vibration 
onset  (right  panel).  C:  Raster  displays  of  discharges.  Each  horizontal  line  corresponds  to  one 
trial  and  each  dot  represents  the  time  of  a  spike’s  occurrence.  Bold  marks  indicate  movement 
onsets  (on  the  left)  and  vibration  onsets  (on  the  right).  The  trials  were  rearranged  in  the  order  of 
increasing  reaction  time  from  top  to  bottom.  Di  Plots  of  the  average  cumulative  sum  (CUSUM). 
The  time  of  statistically  significant  deviations  of  the  CUSUM  trace  are  shown.  E:  Average  wrist 
position  traces.  F:  A  schematic  illustration  of  the  neuron’s  receptive  field  (RF).  This  neuron  had  a 
noncutaneous  RF  deeply  located  in  the  hand.  G:  Cortical  location  of  the  neuron. 

Fig.  6.  Examples  of  characteristics  of  rhythmic  spike  trains  that  seem  unlikely  for  a  renewal 
mechanism  of  rhythm  generation.  A  -  D:  Records  for  a  rhythmically  firing  (~43  Hz)  area  1  neuron 
with  a  multimodal  ISI  distribution.  A:  Activity  during  task  performance  (conventions  as  m  Fig.  5). 
B:  Analyses  of  rhythmic  activity  during  the  hold  period  (conventions  as  in  Fig.  4).  C:  Receptive 
field  (RF)  schematic.  This  neuron  had  a  cutaneous  RF  located  on  the  second  digit.  D:  Cortical 
location  of  the  neuron.  E  -  H:  Records  for  an  area  3b  neuron  that  had  intermpting  spikes.  E; 
Activity  during  task  performance  (conventions  as  in  Fig.  5).  F:  Analyses  of  rhythmic  activity 
during  the  hold  period  (conventions  as  in  Fig.  4).  G:  No  clear  RF  was  found  for  this  neuron.  H; 
Cortical  location  of  the  neuron. 

Fig.  7.  Analyses  of  the  changes  in  activity  of  a  rhythmically  firing  (~30  Hz)  neuron  occurring 
after  vibratory  go-cue  presentation.  A:  A  typical  example  of  vibration-related  changes  in  activity 
(conventions  as  in  Fig.  5).  The  frequency  of  vibration  was  57  Hz.  B:  A  schematic  illustration  of 
the  neuron’s  receptive  field  (RF).  This  neuron  was  activated  by  passive  extension  of  the  fourth 
digit  at  the  metacarpophalangeal  joint.  C;  The  neuron’s  cortical  location.  D:  Analyses  of  the 
rhythmic  activity  during  the  hold  period  (conventions  as  in  Fig.  4).  E:  A  schematic  illustration  of 
the  typical  pattern  of  activity  modulation  by  vibration.  F.  Frequency  distribution  histogram  of 
latencies  of  responses  to  57  Hz  vibratory  stimulation  (the  most  extensively  studied  frequency). 

Fig.  8.  An  example  of  entrainment  of  neuronal  activity  to  the  frequency  of  a  vibratory  stimulus. 

A:  Records  for  a  rhythmically  firing  (~30  Hz)  area  3a  neuron  (conventions  as  in  Fig.  5).  The 
frequency  of  vibration  was  27  Hz.  B:  Phase  raster.  C:  Receptive  field  schematic.  The  neuron 
was  activated  by  passive  wrist  extension.  D:  The  neuron’s  cortical  location.  E:  Analyses  of 
rhythmic  activity  during  the  hold  period  (conventions  as  in  Fig.  4).  F:  Analyses  of  activity  for  the 
epoch  of  sustained  vibratory  response  (100  -  400  ms  after  vibration  onset).  Conventions  as  in  Fig. 
4.  Expectation  and  renewal  density  histograms  were  not  calculated  because  the  duration  of  the 
movement  associated  activity  burst  was  too  short  to  derive  meaningful  displays.  G:  Distribution  of 
discharges  over  the  vibratory  cycle.  The  vibratory  stimulus  is  illustrated  by  the  sinusoidal  trace. 
Upward  deviations  of  this  trace  corresponds  to  upward  movements  of  the  manipulandum. 
Expectation  and  renewal  density  histograms  were  not  calculated  because  the  duration  of  the 
movement  associated  activity  burst  was  too  short  to  derive  meaningful  displays. 
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Fis.  9.  A  typical  example  of  premovement  activity  (PMA.  A:  Records  for  a  rhythmically  firing 
(~44  Hz)  area  2  neuron  (conventions  as  in  Fig.  5).  Vibratory  frequency  was  57  Hz.  B:  No  clear 
receptive  field  was  found  for  this  neuron.  Ci  The  neuron  s  cortical  location.  D.  Analyses  of 
rhythmic  activity  during  the  hold  period  (conventions  as  in  Fig.  4).  E;  A  schematic  illustration  of 
the  typical  PMA  pattern,  i.e.,  an  activity  decrease.  F:  The  frequency  distribution  histogram  of 
PMA  onsets  with  respect  to  movement  onset  and  adjusted  earliest  and  average  EMG  onset  (see 
text). 

Fis  10  Transition  of  neuronal  activity  from  a  rhythmic  to  a  nonrhythmic  pattern  during 
movement.  A:  Records  for  a  rhythmically  firing  (~28  Hz)  area  1  neuron  (conventions  as  m  Fig 
5).  B:  Receptive  field  schematic.  This  neuron  was  activated  by  passive  extension  of  the  second 
digit  at  the  metacarpophalangeal  joint.  C:  The  neuron’s  cortical  location.  D:  Analyses  of  rhythimc 
activity  during  the  hold  period  (conventions  as  in  Fig.  4).  E:  Analyses  of  activity  occurring  with 
movement  (50  -  150  ms  with  respect  to  movement  onset.  Conventions  as  in  Fig.  4.  Full  scale  for 
ISI  histogram  is  80  ms.  Expectation  and  renewal  density  histograms  were  not  calculated  because 
the  duration  of  the  movement  associated  activity  burst  was  too  short  to  derive  meaningful  displays. 

Fig  11  Comparison  of  the  activity  patterns  seen  during  flexion  and  extension  trials.  Records  for 
a  rhythmically  firing  (~38  Hz)  area  3a  neuron  are  shown.  Vibratory  frequency  was  57  Hz  A: 
Records  for  flexion  trials  (conventions  as  in  Fig.  5).  B:  Analyses  of  rhythmic  activity  for  the  hold 
period  during  flexion  trials  (conventions  as  in  Fig.  4).  C:  Schematic  illustration  of  the  neuron  s 
cortical  location.  D:  Records  for  extension  trials.  E:  Analyses  of  rhythmic  activity  for  the  hold 
period  during  extension  trials. 

Fig.  12.  A  model  of  rhythmically  active  SI  neurons  as  tonic  inhibitory  intemeurons  that  gate  the 
activity  of  cortical  output  neurons.  Synapses  marked  black  represent  inhibitory  interactions,  while 
those  marked  white  represent  excitation.  Line  thickness  corresponds  to  the  strength  of  the 
connections. 
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